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We construct the Hubble diagram of Gamma-Ray Bursts (GRBs) with redshifts reaching up to 
z ~ 6, by using five luminosity vs. luminosity indicator relations calibrated with the Cardassian 
cosmology. This model has a major interesting feature: despite of being matter-dominated and 
flat, it can explain the present accelerate expansion of the universe. This is the first study of this 
class of models using high redshift GRBs. We have performed a x-square statistical analysis of 
the GRBs calibrated with the Cardassian model, and also combined them with both the current 
Cosmic Microwave Background (CMB) and Baryonic Acoustic Oscillation (BAO) data. Our results 
show consistency between the current observational data and the model predictions, in particular, 
the best-fit parameters obtained from that x 2 " an alysis are in agreement with those obtained from 
previous investigations. The influence of these best-fit parameters on the redshift at which the 
universe would start to follow the Cardassian expansion, i. e., z ca rd, an d on both the redshift at 
which the universe supposedly had started to accelerate, i. e., z a cc, and the age-redshift relation 
Hoto, are also discussed. Our results also show that the universe, from the point of view of GRBs, 
had undergo a transition to acceleration at a redshift z » 0.2 — 0.7, which agrees with the SNIa 
results. One important point that we notice is that despite the statistical analysis is performed 
with a model that does not need any vaccum energy, we found that the results attained using this 
cosmological model are compatible with those obtained with the concordance cosmology (A-CDM), 
as far as GRBs is concerned. Hence, after confronting the Cardassian scenario against the GRBs 
HD, our main conclusion is that GRBs should indeed be considered a complementary tool to several 
other observational studies for doing precision cosmology. 

PACS numbers: 98.80.-k, 98.70.Rz, 95.36.+X, 98.62.Py 



I. INTRODUCTION 

The discovery of the present late-time acceleration of 
the universe with observations of Type la Supernovae 
(SNIa) [l|, 9, corroborated with the cosmic microwave 
background 0] and the large scale structure observations 
have motivated the introduction of several 
cosmological scenarios. Among them one can quote 
the model with a positive cosmological constant where 
the dark energy evolves with time the model with 
an exotic equation state known as the Generalized 
Chaplygin Gas (GCG), which has the interesting feature 
of allowing the unification of dark energy and dark 
matter [6|, |7|, and similar scenarios. Yet, another 
possible explanation for this accelerate expansion could 
be the infrared modification of gravity as predicted 
by extra-dimension physics, which would lead to a 
modification of the effective Friedmann's equation at late 
times. Such modifications may arise as a consequence 
of our 4-dimensional universe being a surface, or a 
brane, embedded into a higher dimensional bulk space- 
time to which only gravity could spread [1, H, [13, [HI- 
Alternatively, they may arise if there is dark matter self- 
interactions characterized by negative pressure fl2l]. 

One of these interesting possibilities is the modification 
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of the Friedmann's equation by the introduction of an 
additional nonlinear term of mass. Such an idea is 
referred to as the Cardassian model fl3|. 

On the other hand, although SNIa are at the base 
of the suggestion of the late-time acceleration of the 
universe, the fact is that up to the present time it has 
not been possible to register any supernova event at 
redshift z > 2. Therefore, if we wish to know the actual 
expansion history of our universe we need to trace it 
back over a wide range of redshifts. In this perspective, 
the discovery of the X-ray afterglow of the gamma-ray 
burst (GRB) event in 1997 february 28, by the Beppo- 
SAX sattelite [l4[ , allowed the first precise measurement 
of the redshift of a GRB. This breakthrough definitely 
confirmed the long-standing suspicion that GRBs might 
have cosmological origin. That event also makes it 
possible to use GRBs as actual probes in contemporary 
astrophysics and cosmology. This new window onto the 
universe has the advantage of allowing to follow well back 
in time, up to very high redshifts, the expansion history 
of our universe. 

In fact, the possibility of using GRBs as cosmological 
probes has estimulated the search for self-consistent 
methods of bringing them into the realm of cosmology. 
Those procedures include the Amati relation [lfj, the 
Ghirlanda relation [10], the Liang & Zhang relation 
[r^ ]. and the Firmani relation [lg, [13]. All of these 
relations take into account the most relevant physical 
properties of GRBs as the peak energy, jet openning 
angle, time lag and variability. In these lines, recently, 
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Schaefer [20] (hereafter BS2007) used five luminosity vs. 
luminosity indicator relations to calibrate GRBs for a 
specific cosmology. The use of those relations turns 
the GRBs reliable standard candles for practical studies 
in precision cosmology. A similar technique has been 
implemented by Mosquera Cuesta et al. in Ref . [2lj] . 

As is well-known, in the case of supernovae the 
calibration does not depend on any cosmology because 
of a large set of those SNIa are nearby events. In 
practice, however, many of the supernovae in Hamuy et 
al. sample [12] are definitely distant so that the effect of a 
varying cosmology must be introduced, even if it is small. 
Therefore, if we wish to test the various cosmologies up to 
very high redshifts we need to choose a specific cosmology 
and implement the calibration procedure quoted above. 
In this paper we are interested in testing the Cardassian 
cosmology with the current GRBs data after calibrated 
them with the luminosity-distance relation predicted by 
this model. 

This paper is organized as follows: in SectionllTl 
we present an overview of the Cardassian model. In 
Section irm we describe the way the GRBs calibration 
procedure is implemented and demonstrate that it is 
self-consistent by comparing our result for the A-CDM 
scenario with the result obtained by Schaefer (2007). 
Then the construction of the Hubble diagram (HD) of 
those GRBs is made. In Section llVl we perform the best- 
fit analysis for several cases. In Section I VII we study 
other observational aspects of the model like periods 
of both Cardassian and acelerate expansion and the 
age of universe. Finally, in SectionJVlI] we present our 
discussion and conclusions. 



II. 



CARDASSIAN COSMOLOGY 



The dynamics of the universe is worked out through 
both the Friedmann's equation and the evolution 
equations for a perfect, pressureless and non- interacting 
fluid 



H 2 = 

p i + 3H(p i +p i ) = 
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Here the subscript V refers to the components of 
the ideal fluid, a is the scale factor and the subscript 
"0" refers to present-day values. Inspired (apparently) 
on extra dimensions physics, Fresse and Lewis [13 . l23j | 
proposed an explanation for the acceleration phase 
without invoking any vacuum energy or cosmological 
constant, and for a flat universe model as required by 
the CMB observations. Thus, they proposed a modified 
Friedmann equation, Eq.([T]): H 2 = g{p m ), with k = 0, 
where g(p m ) is a different function of the energy density 
and p m contains only radiation and matter (including 



baryon and cold dark matter). The function g(p m ) 
returns to the usual term during the early history of 
the universe, but takes a different form which allows to 
explain the occurrence of an accelerating expansion in 
the recent past of the universe, namely at z ~ 0(1). 
The modified Friedmann's equation reads 



H 2 = Ag(p m ) = Ap n 



Pn 



(l-n) 



(3) 



where A = 87rGjv/3. Thus, for z < z carc i the 
second term dominates [H^. This fact, together with the 
observational data, allows to determine the value of the 
parameter n. In that case, the first term can be neglected 
and thus one finds that the scale factor evolves in time 
following the law: 



a oc £ 3 ™ 



(4) 



so that the (accelerate) expansion is superluminal for 
n < 2/3. The case of n = 2/3 produces a term in the 
Friedmann' equation H 2 cx a~ 2 , which looks similar to a 
curvature term. This feature turns the model attractive 
because the matter alone is sufficient to provide a flat 
geometry. Because of the extra term on the right-hand 
side of the modified Friedmann's equation, the critical 
mass density necessary to have a flat universe can be 
modified so that the total density of the universe (see 
Eqs.(03D) reads 



PtoUz) = PcoidjlOl + zf + Q x f x (z)} . (5) 

Here p to tai{z) = p m (z) + Px(z)- Besides, p c ^i d = 
3Hq/(8itGn) is the usual critical density in the standard 
FLRW cosmology [5!|, and Q°£ s is the observed matter 
density of the universe [60(. The subscript "x" refers 
to any component of the universe that provides an 
additional term in the Einstein's equations. Generically, 
it is called dark energy, but in the Cardassian case it is an 
additional matter term. For the Cardassian model both 
terms in Eq.([5]) come from matter, namely 
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where fij are the fractional density of each component. 
The observed matter density fraction today is given by 
the ratio of the critical mass density of the Cardassian 
universe, p c ,card = Pm,o to that of the standard universe, 
Pc.oid- From Eq.([5]). in the present time we obtain 
directly: 
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" Pc,old [1 + (1 + Zcarf) 3 ^-")] ' 

Conversely, we can express z ca rd, and p C ard in terms of 



a, 



.card = [(n£-- 1 -l] 1/3C1 - B) -l 1 (8) 

oobs r obs-l / n \ 
Pcard = Pc,old" m |" m ~ 1J • (9) 

Finally, from Eq.©, the dimensionless dark energy 
density f x (z) is given by 

/ x (z) = ^J = (l+ Z ) 3 ». (10) 

If the dark energy density corresponds to a 
cosmological constant, one finds that f x (z) = 1, or 
equivalently n — and f2 m = 0.27 for all redshift z. 

Besides, in cosmology the luminosity distance is 
defined as [H,[li] 



where r\ is the co-moving coordinate of the source. 
Using the expression for the propagation of light ds 2 — 
= dt 2 — Yzjri , the modified Friedmann's equation © 
and Eq.© we can re-cast the luminosity distance for a 
flat universe as 

, / x / x c T dz' 

dLiz) = il + z) jr L w)> 

where E 2 {z) = ^ s (l + zf + tt x f x (z). 

For an outlying source of apparent m and absolute 
M magnitudes, distance estimates are made through the 
distance-modulus m—M, which relates to the luminosity 
distance (here expressed in units of Mpc) through: 

fj,(z) =m-M = b\ogd L {z) + 25. (12) 

By plotting the value of the distance modulus, 
fi(z), computed from the estimated luminosity distance, 
dz,(z)), as a function of the redshift (z) one can construct 
the Hubble diagram of gamma-ray bursts after calibrated 
with the Cardassina model. 



III. HUBBLE DIAGRAM OF GAMMA-RAY 
BURSTS 

Gamma-ray bursts are the biggest explosions in the 
universe. The major breakthrough in the understanding 
of GRBs came with the Beppo-SAX satellite discovery of 
the X-ray afterglow of GRB970228 [HI- This allowed the 
first precise determination of the a GRB event redshift, 
what definitely confirmed the long-standing suspicion 
that GRBs had cosmological origin. The huge power 
emitted during a GRB event makes GRBs detectable at 
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FIG. 1: Hubble Diagram of 69 GRBs calibrated with 

the Concordance Cosmology (flat universe with 
Q m = 0.27 and w = -1), using H a = 72 kms _1 Mpc _1 . 
The left panel corresponds to our result (using the 

standard OLS method [2!|) and the right panel 
corresponds to to Schaefer's procedure (Data taken 
from Column (8) of Table 6 in BS2007). 



z~ 20, or even higher, well-deep within the range of the 
epoch of reionization [SIIE HE El| • 

As quoted above, the possibility of using GRBs 
as actual cosmological probes stimulated the search 
for self-consistent methods of bringing GRBs into the 
realm of cosmology. Empirical relations like the Amati 
relation [HI], the Ghirlanda relation [16j |. the Liang & 
Zhang relation [ijj and the Firmani relation [l8| were 
introduced. All of these relations take into account the 
most relevant physical properties of GRBs: the peak 
energy, denoted by -E p0 ak, which is the photon energy 
at which the v F v spectrum is brightest; the jet openning 
angle, denoted by 0j e t, which is the rest-frame time of the 
achromatic break in the light curve of an afterglow; the 
time lag, denoted by r la(5 , which measures the time offset 
between high and low energy GRB photons arriving on 
Earth and the variability, denoted by V, which is the 
measurement of the "spikiness" or "smoothness" of the 
GRB light curve. 

Recently Schaefer in BS2007 presented a subset of 
luminosity vs. luminosity indicator relations that he 
used to calibrate GRBs with A-CDM. He also proposed 
an additional relation between minimum rise time and 
luminosity, where the minimum rise time, denoted by 
t rt , is taken to be the shortest time over which the light 
curve rises by half the peak flux of the pulse. 

For a GRB event to be placed on the HD it is necessary 
to know its isotropic luminosity or its total collimation- 
corrected energy and redshift. The first property can 
not be measured directly but rather it can be obtained 
through the knowledge of either the bolometric peak flux, 
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FIG. 2: Time lag-luminosity relation. The time lags for 
39 GRBs were corrected to the rest frame of each GRB 
event and plotted versus the isotropic luminosity with 
the best-fit (see Eq.([2T])) line superposed on them. 
Notice that the slope in the figure is nearly the 
theoretical value predicted to be —I. 

denoted by Pboio ; or bolometric fluence; denoted by 
SWo (BS2007). Therefore, the isotropic luminosity is 
given by: 

L = 47rd 2 L {z)P holo , (13) 
and the total collimation-corrected energy reads: 

P 7 = 4^|(z)5 bolo F bcam (l + z)- 1 , (14) 

where Fbeam is the beaming factor (1 — cosset). 

The luminosity relations are power-law relations of 
either L or P 7 as a function of r lag , V, £p C ak, t rt . Both 
L and P 7 will be recalculated with luminosity distances, 
and in the case of the Cardassian cosmology with an 
appropriate choice of its cosmological parameters. [6lj 

Table-U below is the same as Table-4 of BS2007. 
Schaefer has compiled all the properties needed to 
calibrate the GRB events, which then produces the GRBs 
HD for any particular cosmology. Column (1) gives the 
six-digit GRB identification number in the usual format 
of YYMMDD. Column (2) gives the redshift (z) of the 
GRB rounded to the nearest 0.01. Column (3) is Pboio 
and its uncertainty in units of erg cm -2 s _1 . Column 
(4) gives 5boio and its uncertainty in units of erg cm~ 2 . 
Column (5) gives the beaming factor, .Fbcam- Column 
(6) gives the time lag for each burst in units of seconds 
in the Earth rest frame and its uncertainty. Column (7) 
gives the variability, V, and its uncertainty. Column (8) 
lists the observed -Speak and its uncertainty, and finally 
Column (9) gives the minimun rise time, t rt , in seconds 
in the Earth rest frame. 



FIG. 3: Variability-luminosity relation. The variability 
for 51 GRBs have been corrected to rest frame of each 
GRB event and plotted vs. the isotropic luminosity 
with the best-fit line (see Eq |2"3")) . This relation was 
verified by Schaefer (35|. 



A. Luminosity Relations and Calibration 
Procedure 

The relationship between a measurable observable of 
the light curve (luminosity indicator) with the GRB 
luminosity is given by the relation luminosity in the form 
of a power-law, i. e., L = B\ && t^ s , L = B v V av ,L — 
B P c^E;ZI\E 7 = B 7 , P oak^ a r k and L = B rt t^ . 
The observed (on Earth) luminosity indicators will have 
different values from those that would be abscribed in the 
rest frame of the GRB event. That is, the light curves 
and spectra seen by the Earth-orbiting satellites suffer 
time dilation and rcdshifting. Therefore, the physical 
connection between the indicators and the luminosity 
in the GRB rest frame must take into account the 
observed indicators and correct them to the rest frame 
of the GRB. For the temporal indicators, the observed 
quantities must be divided by 1 + z to correct the 
time dilation. The observed y-value must be multiplied 
by 1 + z because it varies inversely with time, and 
the observed E pca _k must be multiplied by 1 + z to 
correct the redshift dilation of the spectrum. We have 
also used the same values employed by BS2007 for the 
luminosity indicators to minimize correlations between 
the normalization constant and the exponent during the 
fitting, i. e., for the temporal luminosity we use 0.1s, for 
the variability 0.02, and for the energy indicator 300 keV. 
Notice that these values are only valid for this data. 
Therefore, for another set of GRBs data we can not use 
them. 

To explain the calibration procedure in general, we 
denote the five luminosity relations by Ri — Bi Qf* and 
we take their logarithms to express them as a linear 
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TABLE I: Observables attained from Schaefer's 69 GRBs sample 
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± 0.16 


050126 


1.29 


1.1E-7 


± 


1.3E-8 








2.10±0.30 


0.0039±0.0015 


47 tf 


3.90 


± 0.80 


050318 


1.44 


5.2E-7 


± 


6.3E-8 


3 


46E-6 ± 3.5E-7 


0.0020± 0.0006 




0.0071±0.0009 


47 + 15 

-8 


0.38 


± 0.05 


050319 


3.24 


2.3E-7 


± 


3.6E-8 










0.0028±0.0022 




0.19 


± 0.04 


050401 


2.90 


2.1E-6 


± 


2.2E-7 








0.10±0.06 


0.0135±0.0012 


118±i S s 


0.03 


± 0.01 


050406 


2.44 


4.2E-8 


± 


1.1E-8 








0.64±0.40 




25 tH 


0.50 


± 0.30 


050408 


1.24 


1.1E-6 


± 


2.1E-7 








0.25±0.10 




0.25 


± 0.08 


050416 


0.65 


5.3E-7 


± 


8.5E-8 












1 5 + 23 
iD -2.7 


0.51 


± 0.30 


050502 


3.79 


4.3E-7 


± 


1.2E-7 








0.20±0.20 


0.0221±0.0029 


93 +»» 


0.40 


± 0.20 


050505 


4.27 


3.2E-7 


± 


5.4E-8 


6 


20E-6 ± 8.5E-7 


0.0014± 0.0007 




0.0035±0.0019 


70 tlf 


0.40 


± 0.15 


050525 


0.61 


O.ZHj-D 


_|_ 


7.2E-8 


2 


59E-5 ± 1.3E-6 


0.0025± 0.0010 


0.11±0.02 


0.0135±0.0003 


81 ill 


0.32 


± 0.03 


050603 


2.82 


9.7E-6 


± 


6.0E-7 








0.03±0.03 


0.0163±0.0015 


344tf 2 


0.17 


± 0.02 


050802 


1.71 


5.0E-7 


± 


7.3E-8 










0.0046±0.0053 




0.80 


± 0.20 


050820 


2.61 


3.3E-7 


± 


5.2E-8 








0.70±0.30 




246+™ 

— 40 


2.00 


± 0.50 


050824 


0.83 


9.3E-8 


± 


3.8E-8 












11.0 


± 2.00 


050904 


6.29 


2.5E-7 


± 


3.5E-8 


2 


00E-5 ± 2.0E-6 


0.0097± 0.0024 




0.0023±0.0026 


436+SS 


0.60 


± 0.20 


050908 


3.35 


9.8E-8 


± 


1.5E-8 












41 +_l 


1.50 


± 0.30 


050922 


2.20 


2.0E-6 


± 


7.3E-8 








0.06±0.02 


0.0033±0.0006 


198lf 2 


0.13 


± 0.02 


051022 


0.80 


1.1E-5 


± 


8.7E-7 


3 


40E-4 ± 1.2E-5 


0.0029± 0.0001 




0.0122±0.0004 


olu -20 


0.19 


± 0.04 


051109 


2.35 


7.8E-7 


± 


9.7E-8 












laitJI" 


1.30 


± 0.40 


051111 


1.55 


3.9E-7 


± 


5.8E-8 








1.02±0.10 


0.0024±0.0007 


3.20 


± 1.00 


060108 


2.03 


1.1E-7 


± 


1.1E-7 










0.0032±0.0058 


65 tT ° 


0.40 


± 0.20 



"We take this table of BS2007. The uncertainties given in square 
brackets are conservative estimates for the case in which no error 
bar is quoted in the original literature. 
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TABLE I: Continued ... 

Pbolo Sb l T lag E p<1 ak a TRT 

GRB z (org cm" 2 s" 1 ) (erg cm" 2 ) F b ea m (a) V (keV) (s) 



(1) 


(2) 


(3) 


(4) 


(5 




(6) 


(7) 


(8) 


(9) 




o.oo 


1 TT? 7 -t- 1 fiT7 S 












69 + 19 
02 _ 6 


n An + n 9n 

U.1U It U.iU 


060116 


6.60 


2.0E-7 ± 1.1E-7 












139+t? 

— 36 


1.30 ± 0.50 


060124 


2.30 


1.1E-6 ± 1.2E-7 


3.37E-5 ± 3.4E-6 


0.0021± 


0.0002 


0.08±0.04 


0.0140±0.0020 


237+™ 


0.30 ± 0.10 


060206 


4.05 


4.4E-7 ± 1.9E-8 








0.10±0.10 


0.0025±0.0016 


75 + 12 

10 -12 


1.25 ± 0.25 


060210 


3.91 


5.5E-7 ± 2.2E-8 


1.94E-5 ± 1.2E-6 


0.0005± 


0.0001 


0.13±0.08 


0.0019±0.0004 


149lJ°° 


0.50 ± 0.20 


060223 


4 41 


2.1E-7 ± 3.7E-8 








0.38±0.10 


0.0075±0.0033 


?1 +100 


0.50 ± 0.10 


060418 
060502 


1.49 
1.51 


1.5E-6 ± 5.9E-8 
3.7E-7 ± 1.6E-7 








0.26±0.06 
3.50±0.50 


0.0070±0.0005 
0.0010±0.0017 


230+gj 


0.32 ± 0.08 
3.10 ± 0.30 


060510 


4.90 


1.0E-7 ± 1.7E-8 










0.0028±0.0019 


n= +[60] 
M0 -[30] 




060526 


3.21 


2.4E-7 ± 3.3E-8 


1.17E-6 ± 1.7E-7 


0.0034± 


0.0014 


0.13±0.03 


0.0112±0.0039 


25 +[5] 


0.20 ± 0.05 


060604 


2.68 


9.0E-8 ± 1.6E-8 








5.00±1.00 




40 +1 


0.60 ± 0.20 


060605 


3.80 


1.2E-7 ± 5.5E-8 








5.00±3.00 




16 „+[200] 


2.00 ± 0.50 


060607 


3.08 


2.7E-7 ± 8.1E-8 








2.00±0.50 


0.0059±0.0014 


120±if 


2.00 ± 0.20 



"We take this table from BS2007. The uncertainties given in 
square brackets are conservative estimates for the case in which no 
error bar is quoted in the original literature. 



relation of the form 

log Ri = log Bi + a. L log Q t y l = b + a x. L . 

In order to use the linear regression method to 
determine the best-fit parameters we also take into 
account that in those luminosity relations these variables 
are independent. This feature allows to use the OLS 
Bisector method [2{|- Besides, because the scatter of the 
data is consistent with a Gaussian distribution we can 
use the error propagation law: The optimum f for the 
quantity of interest /(ctj) i = 1, 2, ... n calculated from 
the medians di and their standard deviation aj is given 




Then, for the luminosity indicators, applying this law, 
we obtain the standard deviation associated to the x- 
variable 




Then, the calibration will essentially be a linear fit on 
a log-log plot of the luminosity indicator Qi versus the 
burst luminosity Ri, i. e., we fit the line yi = b + axi, 
and find the corrected value of the luminosity using 
the subroutine Sixlin.f [liH upon computing y~i(a,b,Xi). 
Therefore, their associate error or standard deviation, 
using the Eq. p5)l . is given as 

4 =<?l + {<?bXi) 2 + {b<J Xl ? , (18) 

where o~ Xi is given by Eq. (|16[) . As the uncertainties in 
both luminosities and their indicators are small than 
the observed scatter about the best-fit line, we must 
introduce an additional source of intrinsic scatter in the 
last expression, denoted by cr Bya . This value can be 
estimated by finding the value such that a x 2 fit of the 
luminosity calibration produces a reduced \ 2 of unity 



dx(Qi 
dQi 



'Qi 



1 Q-Qj 
In 10 Q l 



Similarly, the standard deviation associated with the 
y-variable is function of -Pbolo when we use the burst 
luminosity Eq.fT3]), or is function of S'boio and Fbeam 
if we use the total collimation-corrected energy Eq.dll 

In the first case 



1 (T P 



In 10 Pbolo 



(16) 



For the second case, corresponding to log E 1 , the 
standard deviation is given by 



X 



N (- 



i 



N 4- o& + <? 2 



(19) 



Here cr|. is given by Eq. (fl"8]) and then we re-cast the 
standard deviation for the correct luminosity as: 



+ {<T b x l ) 2 + (ba Xi ) 2 + o~l 



(20) 



With the purpose of verifying our procedure we used 
the standard OLS method [29| in a barely different form 
to that one used by Schaefer (see pag.26 of BS2007), 
and then we calibrated the GRBs for the Concordance 
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FIG. 4: iJpoak-luminosity relation. The -E pca k values for 
64 GRBs have been corrected to the rest frame of each 
GRB event and plotted versus the isotropic luminosity 
with the best-fit line, (see Eq l2"5"l) . This relation was 
proposed and verified by Schaefer [20l |40| . 



Cosmology (flat universe with O m = 0.27 and w = —1). 
We choose H = 72 kms^Mpc -1 from the HST Key 
project [32]. In FigfTJwe show our results and we compare 
them with the HD obtained by Schaefer (2007). One 
can notice that the difference between both plots is very 
small, what confirms that our method is self-consistent. 

In the next subsections we present the results found 
with the calibration procedure of the five luminosity 
relations. All of them are based on an assumed 
Cardassian model using Eq. (Tr2|) with f2 m = 0.27, n = 0.2 
and H = 72 kms^Mpc- 1 !^ 



value of logL using Eq. (|2"0]) (the of agsys -value found is 
0.36) is given by 



°fogL = a l + i a b 



(l + z)0.1s 



0.4343 b criag 



' "fag.sys • ( 22 ) 



2. Variability versus Luminosity 

The variability (V) of a GRB event was identified as 
indicator luminosity by Fenimore & Ramirez- Ruiz [36l ]. 
Subsequently, Reitchart [371 ] proposed a new relation 
between variability and isotropic luminosity, which is 
similar to the time lag-luminosity. Its origin is based 
in the physics of the relativistic shocked jets[H, l39j . 
The calibration plot for the V — L relation is given in 
FigEl along with the best-fit line. This best-fit can be 
represented by the equation 



logL = 52.43 + 1.77 log 



V(l + z) 



0.02 



(23) 



The l-cr uncertainties in the intercept and slope are 
a a = 0.07 and <Jb — 0.19, the ov, sys value found was 0.47 
and the uncertainty in the log of the luminosity, using 
Eq.® is 



2 2 
^logL = a a- 



V(l + z) 



0.02 



0.4343 6ct, 
V 



v,sys 

(24) 



Epeak versus Luminosity 



Time 



versus Luminosity 



The time lag, ri ag , was identified as a luminosity 
indicator by Norris et al. [33| who proposed a power 
law relation between the time lag and the isotropic 
luminosity. This relation is a consequence of the 
empirical/theoretical Liang-Kargatis relation [34| and 
was verified by Schaefer [35|. The calibration data are 
plotted in FigOas logri ag /(l + z) versus L. The best-fit 
linear regression line is also plotted. The equation for 
this calibration line reads: 



logi = 52.23- 1. 00 log 



(1 + z)0.1s 



(21) 



Notice that the slope value satisfactorily agrees with 
the theoretical value of -1. The l-a uncertainties in 
the intercept (a = 52.23) and slopes (b = —1.00) are 
a a = 0.07 and o\> — 0.09, and the uncertainty on the 



This luminosity relation was proposed by Schaefer [401] 
and is related to the instantaneus physics at the time of 
the peak. It was also verified in BS2007. The calibration 
plot for the E vea k — L relation is given in FigfJ] along with 
the best-fit line. This best-fit line can be represented by 
the equation 



logi = 52.18+ 1.68 log 



-Epeak (1 + Z) 



300 keV 



(25) 



The l-cr uncertainties in the intercept and slope are 
a a = 0.05 and erf, = 0.10, and the <T_B pcak , S ys value is 0.4. 
The uncertainty in the log of the luminosity using Eq. (|20[) 



2 



Ob 



E, 



peak 



(1 + *) 



300 keV 



0.4343 6 (tr 



E, 



peak 



pcak,sys 



(26) 




FIG. 5: -E pea k-E 7 relation, or Ghirlanda relation. The 
-Epeak values for 27 GRBs have been corrected to the 
rest frame of the GRB event and plotted versus the 
total burst energy in the 7-rays with the best-fit line 
(see Eql27|). This relation allows an empirical correction 
to the determined luminosity distance for each GRB 
event, in the same way as the light-curve shape 
correction for SNIa. 



FIG. 6: Minimum rise time-luminosity relation. The 
rise time for 62 GRBs have been corrected to the rest 
frame ofthe GRB event and plotted versus the isotropic 
luminosity with the best-fit line (see Eq|2"9"|). This 
relation was proposed and confirmed for Schaefer in 
BS2007. 

5. Rise Time versus Luminosity 



4- Epaak versus E 1 

Ghirlanda et al. [l(| found that for GRBs the total 
energy emitted in 7-rays (E 7 ) after a proper collimated 
correction, correlates tightly with the peak energy E'peak 
(in the v-¥ v spectrum). Thus the isotropically equivalent 
burst energy could be determined with sufficient accuracy 
to be used in a practically fashion for cosmological 
studies. The physics of this relation is explained within 
the standard jet model [IE [H, [H, EH- The calibration 
plot for the -Epeak-E-y relation is given in Fig|5] along with 
the best-fit line. This best-fit can be represented with the 
equation 



\ogE 1 = 50.52 + 1.68 log 



Epeak (1 + Z) 

300 keV 



(27) 



The 1-er uncertainties in the intercept and slope are 
u a = 0.05 and a b = 0.10, and 0£ poakiSys = 0.21. The 
uncertainty in the log of the luminosity is obtained using 
the Eq.® is 



In BS2007, in an effort to understand the physical 
origin of the variability, Schaefer calculated the 
variability for a wide range of simulate light curves 
constructed from individual pulses. He found that the 
most important determinant of the V-value was the 
rise time in the light curves, and this rise time can be 
connected with the physics of the shocked jet. The 
calibration plot for the r RT — L relation is given in 
FigJS] along with the best-fit line. This best-fit can be 
represented with the equation 



logL = 52.48- 1.21 log 



(1 + z) 0.1s 



(29) 



The l-cr uncertainties in the intercept and slope are 



<j a = 0.07 and cjb = 0.11, and the value of cxrt.; 



0.47. 



The uncertainty in the log of the luminosity is given by 
Eq.© 



7 logL 



Ob 



(l + z)0.1s 



2 



-Epeak (1+3) 



300 keV 



0.4343 b g e 



E, 



peak 



ipcak.sys 



0.4343 b a T 



(30) 



In Table-fll] we collect our results. Notice that they 
(28) are very similar to the results obtained in BS2007 for 
the Concordance Cosmology (57 m = 0.27 and w = 



TABLE II: Calibration Results 



Luminosity Relation 


a 




b 


<?b 


cr ays 


"Hag — L 


52.23 


0.07 


-1.00 


0.09 


0.36 


V - L 


52.43 


0.07 


1.77 


0.19 


0.47 




52.18 


0.05 


1.68 


0.10 


0.40 


^peak~E~y 


50.52 


0.05 


1.62 


0.10 


0.21 


7rt — L 


52.48 


0.07 


-1.21 


0.11 


0.47 



— 1) and for the dark energy parametrization (wo = 

— 1.31,to =1.48). Indeed, the set of x 2 for both methods 
are as follows: x 2 A - CDM/Schafer = 71.61248, 
X 2 A - CDM/Schafer dof = 1.03786, x 2 A - CDM/Ours = 
72.16559, x 2 A - CDM/Ours dof = 1.04587, x 2 Card = 
72.034, x 2 Carddof = 1.04379. This confirms one more 
time that the calibration procedure depends weekly of 
the input cosmology. 

6. Combining the Derived Distance Moduli 

After performing the calibration procedure for every 
L, or Ey, we can also derive a luminosity distance, 
or equivalently their logarithms, and estimate their 
associate distance moduli using Eqs. (|12ll3ll4[) obtaining: 

M (z) = ^log(L)-^log(4^P bolo ) + 25, (31) 

n{z) = - log (b 7 ) - - log I — y^~ z — ) • ( 32 ) 

The propagate uncertainties will depend on whether we 
calculate the distance modulus using L or Ej. Therefore, 
using Eq. (|15p . one can easily obtain 

<r M = 2.5 [a 2 + a 2 ] 1/2 , (33) 

where a y takes either the form of Eq. (fTo| if we calculate 
L, or the form of the Eq. (fl"7)) if we calculate E 1 , and ay 
is given by Eq. (|20[) applied to the respective luminosity 
relation. 

With five luminosity indicators, each burst will have 
up to five estimated distance moduli an their 1-er 
uncertainties. As in BS2007, we label as [i\ ± , 
(12 ± cr M2 , ^ 3 ± cr M3 , ± tT P4 , /i 5 ± 0" M5 , for the five 
indicators, respectively. Then to calculate the best 
estimate fi for each GRB event it is necessary to make a 
weight-average of all available distance moduli, which is 
given by: 



H = — WiHi , (34) 

i 

where Wi = 1/of and w = w i 




Redshift (z) 



FIG. 7: Hubble Diagram of 69 GRBs after having been 
calibrated with Cardassian Cosmology. The blue line 
corresponds to the HD predicted by the Cardassian 
parameters used in this calibration; f2 m = 0.27, n = 0.2. 
The black line corresponds to the HD for the 
Concordance Cosmology (A-CDM model with 
f2 m = 0.27, w = —1). Here we assumed 
H = 72kms- 1 Mpc- 1 . 

Then applying the error propagation law (see Eq |15p 
to Eq. (|34|) we obtain the standard deviation associate to 
this best estimate as <r M = 1/^/w. 

B. Hubble Diagram 

In Table- lnil we collect the distance moduli calculated 
for each burst. Columns (1) and (2) give the GRB event 
and their redshift (z), Column (3)-(7) give the distance 
modulus fii ± <7 Mi calculated for each luminosity relation 
for the Cardassian cosmology: £l m = 0.27, n = 0.2 
and Ha = 72 kms _1 Mpc _1 , and the Column (8) gives 
the combined distance modulus values obtained with the 
Eq. (j34j) and their uncertainties. 

Fig[7] shows the HD for the GRBs calibrated with the 
Cardassian Cosmology. Notice that the observational 
data points are consistent with the Cardassian model 
predictions because its %-square per degree-of-freedom 
X-dof ~ 1-00 ~ 1-07 (see next section). 



IV. DATA ANALYSIS: THE METHOD 

We can now determine the best-fit to the set of 
cosmological parameters 8 = (Q nl ,n) after imposing 
several constraints. Later on we combine GRBs data with 
those from the CMB and BAO to put tigh constraints 
on these parameters. To this end, we construct the log- 
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TABLE III: Derived Distance Moduli 







Mi 


A*2 


^3 


H4 


Ms 


A 1 


GRB 


z 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


970228 


0.70 




42.41±1.23 


42.31±1.19 




43.29±1.20 


42.67±0.70 


970508 


0.84 


43.00±1.14 


42.99±1.26 


45.54±1.05 


43.78±0.66 


42.94±1.21 


43.77±0.43 


970828 


0.96 




42.85±1.21 


43.96±1.05 


43.47±0.66 


43.14±1.24 


43.43±0.47 


971214 


3.42 


48.63±1.46 


48.54±1.22 


47.44±1.04 




49.17±1.32 


48.31±0.62 


980613 


1.10 






45.75±1.35 






45.75±1.35 


980703 


0.97 


44.42±0.97 


44.81±1.20 


45.98±1.04 


43.45±0.61 


42.00±1.25 


44.02±0.41 


990123 


1.61 


43.14±0.95 


44.69±1.20 


45.49±1.05 


45.38±0.69 




44.78±0.46 


990506 


1.31 


44.69±1.09 


44.08±1.19 


44.07±1.04 




43.80±1.21 


44.18±0.56 


990510 


1.62 


46.45±1.02 


45.11±1.19 


44.13±1.03 


45.38±0.60 


45.53±1.20 


45.34±0.41 


990705 


0.84 




45.11±1.20 


43.47±1.03 


42.85±0.66 


45.66±1.30 


43.69±0.47 


990712 


0.43 






41.75±1.07 


41.41±0.69 




41.51±0.58 


991208 


0.71 




40.38±1.22 


42.09±1.04 




41.88±1.20 


41.52±0.66 


991216 


1.02 


43.43±1.01 


42.38±1.19 


42.61±1.04 


43.53±0.68 


43.19±1.23 


43.16±0.43 


000131 


4.50 




46.96±1.22 


47.59±1.04 




48.34±1.37 


47.58±0.69 


000210 


0.85 




40.78±1.22 


43.68±1.03 




41.81±1.21 


42.28±0.66 


000911 


1.06 




44.39±1.21 


45.54±1.06 




44.66±1.31 


44.94±0.68 


000926 


2.07 




46.12±1.22 


44.13±1.03 




47.23±1.44 


45.49±0.69 


010222 


1.48 




43.20±1.19 


43.56±1.03 


44.90±0.57 


43.55±1.23 


44.28±0.43 


010921 


0.45 


42.76±1.02 


40.86±2.42 


43.07±1.07 




41.06±1.26 


42.34±0.62 


012111 


2.14 






46.96±1.06 


44.98±0.67 




45.55±0.57 


020124 


3.20 


47.73±1.18 


48.37±1.28 


46.14±1.07 


46.39±0.67 


47.22±1.23 


46.88±0.44 


020405 


0.70 




43.90±1.19 


44.40±1.12 


43.41±0.79 


42.56±1.21 


43.55±0.52 


020813 


1.25 


44.31±0.97 


45.19±1.19 


43.91±1.04 


43.90±0.64 


42.86±1.21 


44. 01±0. 41 


020903 


0.25 






40.65±1.29 






40.65±1.29 


021004 


2.32 


46.34±1.21 


46.60±2.76 


46.62±1.18 


45.71±0.84 


47.53±1.34 


46.34±0.53 


021211 


1.01 


43.98±0.94 




42.19±1.06 




44.45±1.21 


43.51±0.61 


030115 


2.50 


46.48±1.09 


47.25±1.78 


46.42±1.14 




45.36±1.29 


46.29±0.63 


030226 


1.98 


46.84±1.44 


47.07±1.97 


46.64±1.09 


46.10±0.69 


46.35±1.26 


46.38±0.48 


030323 


3.37 






46.74±1.58 




47.22±1.46 


47.00±1.07 


030328 


1.52 


45.13±1.43 


44.61±1.22 


44.84±1.04 


44.49±0.64 




44.65±0.47 


030329 


0.17 


41.94±0.98 


41.55±1.20 


39.57±1.03 


38.81±0.61 


40.49±1.21 


39.97±0.41 


030429 


2.66 




47.65±2.33 


45.44±1.14 


46.44±0.89 


46.57±1.26 


46.28±0.59 


030528 


0.78 


42.75±1.05 


44.75±2.08 


44.20±1.10 




46.08±1.26 


44.19±0.62 


040924 


0.86 


43.83±0.95 




42.60±1.04 




45.09±1.20 


43.74±0.61 


041006'' 


0.71 




44.09±1.19 


42.38±1.10 


44.08±0.73 


43.26±1.24 


43.60±0.50 


050126 


1.29 


45.37±0.98 


46.74±1.41 


45.77±1.08 




44.69±1.27 


45.57±0.57 


050318 


1.44 




46.33±1.22 


44.20±1.08 


45.78±0.72 


46.14±1.21 


45.60±0.49 


050319 


3.24 




46.49±1.93 






48.65±1.23 


48.03±1.04 


050401 


2.90 


46.06±1.15 


46.94±1.22 


45.22±1.06 




48.56±1.30 


46.52±0.59 


050406 


2.44 


48.16±1.19 




46.41±1.43 




48.96±1.45 


47.87±0.77 


050408 


1.24 


45.17±1.05 








45.76±1.27 


45.40±0.81 


050416 


0.65 






41.38±1.12 




45.22±1.43 


42.83±0.88 


050502 


3.79 


47.26±1.46 


50.00±1.30 


46.89±1.27 




47.16±1.39 


47.87±0.67 


050505 


4.27 




46.97±1.59 


46.86±1.21 


48.43±0.99 


47.60±1.30 


47.63±0.61 


050525 


0.61 


44.01±0.95 


44.26±1.19 


41.93±1.03 


43.13±0.71 


43.32±1.19 


43.27±0.43 


050603 


2.82 


45.69±1.45 


45.60±1.23 


45.48±1.07 




44.60±1.20 


45.32±0.61 


050802 


1.71 




45.74±2.52 






45.34±1.25 


45.42±1.12 


050820 


2.61 


45.87±1.05 




48.43±1.08 




44.97±1.26 


46.55±0.65 


050824 


0.83 










43.22±1.37 


43.22±1.37 


050904 


6.29 




47.06±2.48 


51.06±1.13 


49.19±0.77 


47.77±1.27 


49.27±0.55 


050908 


3.35 






46.82±1.06 




46.91±1.23 


46.86±0.80 


050922 


2.20 


46.46±1.01 


43.91±1.25 


45.86±1.04 




46.43±1.21 


45.77±0.56 


051022 


0.80 




43.47±1.19 


44.69±1.03 


43.73±0.57 


43.33±1.22 


43.81±0.43 


051109 


2.35 






46.59±1.11 




44.50±1.27 


45.69±0.83 


050922 


1.55 


44.90±0.96 


44.64±1.35 






43.71±1.30 


44.52±0.67 


060108 


2.03 




46.90±3.83 


46.87±1.52 




48.04±1.74 


47.34±1.09 


060115 


3.53 






47.34±1.04 




48.38±1.36 


47.73±0.83 


060116 


6.60 






49.29±1.28 




47.05±1.42 


48.29±0.95 


060124 


2.30 


46.83±1.09 


47.39±1.24 


46.89±1.10 


46.96±0.69 


46.03±1.27 


46.87±0.45 




41. UD 


A c fi/l -1- 1 A A 
± . 414: 




af. nfi 

410.OO3I 1 .UO 




A K 71-1-1 99 


AF. c;Q_i_n cc 
410. Do IU. DO 


060210 


3.91 


47.48±1.15 


45.08±1.27 


47.52±1.11 


49.43±0.73 


46.63±1.30 


47.84±0.46 


060223 


4.41 


47.47±0.99 


48.94±1.48 


47.39±1.07 




47.80±1.23 


47.74±0.58 


060418 


1.49 


44.90±0.96 


45.19±1.20 


45.99±1.04 




45.24±1.23 


45.33±0.54 


060502 


1.51 


43.60±1.08 


42.99±3.53 


46.81±1.19 




43.79±1.31 


44.65±0.67 


060510 


4.90 




48.02±1.77 


48.89±1.19 






48.62±0.99 


060526 


3.21 


48.22±0.98 


49.09±1.38 


44.88±1.10 


46.83±0.82 


48.53±1.24 


47.29±0.47 


060604 


2.68 


45.16±1.01 




46.56±1.07 




47.98±1.28 


46.36±0.64 


060605 


3.80 


45. Mil. 26 




49.37±1.19 




46.44±1.34 


47.11±0.73 


060607 


3.08 


45.08±1.03 


47.67±1.31 


47.56±1.10 




45.35±1.25 


46.33±0.58 



b Notice that in BS2007 this GRB has not associated Ti ag , but he 
obtained a /ii-value for this GRB event. 
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TABLE IV: Best-fit analysis 



Case c £l m ±aa m d n ±a„ d \ 2 xlof 



exact H 0.45 ±0.20 [±o.J|] -0.70 ±2.68 [t 00 ] 70.88 1.06 

exact H , prior fi m 0.28 ±0.04 [to.ot] 0.27 ±0.16 [tj^g] 71.53 1.07 

Prior Ho 0.40 ±0.36 [to.H] 0.26 ±1.66 [ttool 73 -28 1.09 

Prior H +TZ 0.35 ±0.09 [±S:os] °- 39 ± - 22 toils] 73 -29 1.09 

Prior Ho.+A 0.33 ±0.04 [1°;°*] 0.45 ±0.29 73.30 1.09 

Prior # .±ft±.4 0.32 ±0.04 0.34 ±0.16 [tg^f] 73.50 1.10 



c Here we refer to prior a free parameter with Gaussian 
Prior distribution. The marginalization procedure takes Hubble 
parameter Ho as nuissance parameter. 

d The parameter errors given in square brackets are estimates 
obtained by the MINOS analysis of the base code MINUET, 
corresponding to one-standard deviation. En the case that they 
are zero, the MENOS can not determine such errors. 



likelihood function x 2 by assuming hypothetical m and 
M magnitudes fitted to a number of redshifts: 



terms of the ffrj-hidependent luminosity distance: Dl 
HqcIl (where cLl is the proper luminosity distance) [46| 



Xo(0,#oA) 



^ [Vte{Zi,0,H O ) 
i=l 



where N is the number of data, Hd.a.t{z) and cr^ dat (^) 
are, respectively, the fitted distance modulus and its 
dispersion at redshift z obtained from the calibration 
procedure. In addition, fith(z) is the theoretical 
prediction of the distance modulus and bi(i = 1, ... ,5) 
are the slopes obtained from the regression analysis. 

We can obtain the function x 2 in the parameter 
space (9, Hq) by marginalizing over the five nuisance 
parameters 6j. However, a reasonable approximation is 
to consider the values of the slopes fixed at the values 
obtained in the regression analysis, because the GRB 
HD is nearly independent of the assumed cosmology. 
Therefore, the x 2 function reads 



2/7) it \ (pth(zi,0,H ) - nnt(zi)) 2 

Xo(V, Hq) = ^ 



Mi) 



(36) 



where /xg t and <J^ iit are obtained in the previous section 
(see Eq.flM]) and Column (8) in Table-|TTI|- 

On the other hand, if we want to treat Q m as a free 
parameter with Gaussian prior distribution centered in 



Qtrue w jt n spread an n 



-prior 



then we have to use 



truc\2 



(37) 



-prior 



A. GRBs bounds ± constraints from other 
cosmological observations 

To combine the GRBs data with the CMB we consider 
the model-independent shift-parameter 1Z, defined in 



K 



1 + Z T 



dz 

W) 



(38) 



The shift-parameter 1Z is related to the position of 
the first acoustic peak of the CMB anisotropy power 
spectrum but it can not be directly measured. However, 
its value is estimated from the data assuming a flat 
cosmology with dark matter and cosmological constant. 
The shift-parameter is defined, approximately, as the 
ratio of the sound horizon at recombination to the co- 
moving distance to the last scattering. The quantity 
E(z) is given by Eq. ([T^|) and z r = 1089 is the redshift of 
recombination. (From the three-year results of WMAP 
[47l | Wang and Mukherjee estimated 1Z= 1.70 with spread 
(t-ji = 0.03 [48]). Then, the statistical significance of a 
model is determined evaluating x\ = — l-70) 2 /cr^ 
together with the x 2 01 the GRBs data, i. e., 



x 2 = xl{o,H,) + x 2 n {e) 



(39) 



Meanwhile, to combine the GRBs data with the BAO 
we considered the ^-parameter which describes the BAO 
peak in the matter perturbation power spectrum, given 
by 0: 



A = vn 



dz' 



2/3 



(40) 



where z\ = 0.35 is the redshift at which the acoustic 
scale has been found. Eisenstein [5§] found that the 
estimate value of the ,4-parameter is A— 0.469 with 
spread 04 = 0.017. Similarly to the shift-parameter, 
the statistical significance of a model is then determined 
evaluating Xa — (A — 0.469) 2 /cr^ together with the x 2 
of the GRBs data, i. e., 



x 2 = x 2 o(o,H ) + x 2 A (e) ■ (41) 

This BAO peaks are clearly seen in the matter power 
spectrum, which together with the overall shape put tight 
constraints on the model parameters. 

Therefore, if we combined the GRBs data plus CMB 
and BAO data the most general expression to x 2 is given 
by: 

X 2 (9,H ) = x 2 o(0,H ) + x 2 z(0) + X 2 a(0) ■ (42) 

To get the Q m — n parameter space we use the 
marginalization procedure choosing Ho with Gaussian 
Prior distribution as nuissance parameter. Therefore, we 
can define a modified x 2 -statistic by: 




X 2 (9) = -2\n 



where tt{Hq) 



dH exp[-0.5x 2 (e,H )}ir(H ) 



(43) 



Here the 



function x 2 (9,H ) takes the form of Eqs.©, ©, SB) 
or (|42|) depending on the analysis that will be performed. 
It is straightforward to minimize x 2 (6*) using the base 
code MINUIT[50[ or the function FindMinimum of the 
Software Mathematical to find Xmin(^)> where Xmin * s 
the minimum obtained for the best-fit parameters values 
§=(Q m ,n). 

The variable Xmin ^ s random in the sense that it 
depends on the random dataset that is used. Its 
probability distribution is a x 2 distribution for N — v 
degrees of freedom (y is the number of free parameters). 
This implies that the 68% of the random parameters in 
a given dataset will give a x 2 such that 



FIG. 8: Confidence- level contours in the VL m — n 
parameter space for the flat Cardassian Model. The 
solid lines separate the 1-er, 2-cr and 3-cr confidence-level 
(color) regions for the sole GRBs data, where the black 
point (0.45, -0.70) corresponds to the best fit. The 
encircling dashed lines divide the l-cr, 2-cr and 3-cr 
confidence-level regions obtained after imposing the 
prior knowledge on fi m , with the green point (0.28, 
0.27) corresponding to the best-fit value, while the blue 
point (0.27, 0) corresponds to the cosmological constant 
A. One can see that in both cases the cosmological 
constant lies in the 2-cr level. 



cosmological model under consideration. We divide the 
statistical analysis in two subsections, considering the 
GRBs data alone and then combining GRBs with CMB 
and BAO data. 



X 2 -xL„< Axi», (44) 

where Ax 2 a (^) is 2.3 for two [y — 2) free parameters. 
Thus, Eq. pi)) defines the 1-cr surface around the best-fit 
parameter. Similarly, it can be shown that 95.4% of the 
random numbers in the dataset will give a x 2 such that 

X 2 -xL„<AxLH, (45) 

where Ax2 (T ( l/ ) is 6.17 for two free parameters. Thus, 
Eq. pS")) defines the 2-cr surface around the best-fit and 
similarly for higher cr's. 

The fixed values of Ax 2 CT (i = 1,2, . . .) determine the 
boundary of the i-th confidence-level contours and gives 
the probability of having a given x 2 for the true values of 
the estimate parameters laying inside the boundary. This 
criterion permit us to reject or accept any pair of (f2 m , n) 
parameter, and if X 2 /{N— v) w 1 the fit is said to be good 
and the data are said to be consistent with the considered 



V. STATISTICAL ANALYSIS: RESULTS 

A. GRBs data alone 

In this subsection we consider only the GRBs dataset 
for the following three cases, and perform their statistical 
analysis and construct their confidence-level contours, 
according to the method sketched above. 

exact Ho 

Here we suppose an exact knowledge of the Hubble 
parameter today. For the minimization procedure we 
use Eq.® obtaining Xmin = 70.88 (Xd„f= L06 )- In 
this case, we have two free parameters, the best-fit value 
corresponds to the point (0.45, —0.7) in the fi ro — n 
parameter space. The confidence contours are shown in 
the FiglH] Hence, we see that the case n = 0, which 




FIG. 9: Idem to FigHJ but now the black point (0.40, 
0.26) corresponds the best-fit value, and the blue point 
(0.27, 0) corresponds to cosmological constant A. One 
can see that in this case the cosmological constant lies 
in the 1-cr level. 



corresponds to the cosmological constant A, lies in the 
2-cr level. 



exact Ho, prior £2 m 

In this case we consider an exact knowledge of Hubble 
constant today and the matter density as free parameter 
with a uniform prior in the range fl m = 0.27 ± 0.04. The 
best- fit analysis was performed using Eq. (|37[) obtaining 
Xmin = 71.53 (xiof = 1-07), and the values correspond to 
the point (0.28, 0.27) in the f2 m — n parameter space. The 
confidence-level contours correspond to the filled regions 
limited by dashed lines in Figj8] In this case, the point 
that corresponds to A also lies in the 2-cr level. 



Prior Ho 

Here the best-fit analysis is the most realistic 
constraint for GRBs. The minimization procedure was 
performed using the Eqs. (|36|) and (|43|) and the Xmin = 
73.28, (Xdof— 1-09). The best-fit value corresponds to 
the point (0.40, 0.26) in the f2 m — n parmeter space. The 
result in the f2 m -value agrees with the first one estimated 
by Schaefer in BS2007. The confidence-level contours 
correspond to the filled regions limited by solid lines 
in FigsOU [TUl HU and [T3J Hence, one can see that the 
cosmological constant lies in the Tcr level. 




FIG. 10: Confidence-level contours in the f2 m — n 
parameter space for the flat Cardassian Model with 
GRBs + TZ. The solid lines separate the 1-cr, 2-cr and 
3-cr confidence-level (color) regions for the GRBs data 
alone, with best-fit value (0.40,0.26) represented by the 
cyan point. The encircling dashed lines divide the 1-cr, 
2-cr and 3-cr color-filled confidence regions of the bounds 
from GRBs + TZ, with the best-fit value (0.35,0.39) 
represented by the green point. One can see that 
cosmological constant A (blue point) lies in the 2-cr 
level. 



B. Combined Data 

In this subsection we analyze the role of the GRBs with 
combined constraints from the CMB and BAO data. 



Prior Ho and TZ 

Here we perform the minimization procedure using 
Eq. (|4*3")) without Xprior an d Xa> an d found Xmin = 73.29, 
(Xdof= 1-09). The best- fit value corresponds to the point 
(0.35, 0.39), and the confidence contours are shown in the 
FigHD] One can see that the relevant constraint on the 
parameter n, that would represent the dark energy and 
the cosmological constant A, lies in the 2-cr level. 



Prior Ho and A 

Here we perform the minimization procedure using 
Eq. (|4*3"l) without Xprior an d Xtj> anc ^ found Xmin = 73.30, 
(Xdof = 1-09). The best- fit value corresponds to the point 
(0.33,0.45) in the f2 m — n parameter space, and the 
confidence contours are shown in Fig[TTJ We can see 
that the cosmological constant A lies in the 1-cr level. 



U.l 0.2 U.j U.4 U.5 0.6 u. / u.s o.y 1 1.1 



U.l 0.2 0.3 U.4 0.5 U.b U. / U.S U.y 1 1.1 1.2 



FIG. 11: Confidence- level contours in the f2 m — n 
parameter space for the flat Cardassian Model with 
GRB plus A. The solid lines separate the 1-cr, 2-cr and 
3-er confidence-level (color) regions for the GRBs data 
alone, with best-fit value (0.40,0.26) represented by the 
cyan point. The encircling dashed lines divide the 1-er, 
2-cr and 3-cr color-filled confidence regions of the bounds 
from GRBs plus A, with the best-fit value (0.33, 0.45) 
represented by the green point. One can see that 
cosmological constant A (blue point) lies in the 1-cr 
level. 



FIG. 12: Confidence- level contours in the Q m -n 
parameter space for the flat Cardassian Model with 
GRBs + 1Z + A. The solid lines separate the 1-er, 2-cr 
and 3-cr confidence-level (color) regions for the GRBs 
data alone, with best-fit value (0.40,0.26) represented 
by the cyan point. The encircling dashed lines divide 
the 1-cr, 2-cr and 3-cr color-filled confidence regions of 
the bounds from the GRBs + 1Z + A, with the best-fit 
value (0.32, 0.34) represented by the green point. Hence, 
one can see that cosmological constant A (blue point) 
lies in the 2-cr level. 



Prior Ho, 1Z and A 



Here we consider the most general case for constraining 
the cosmological parameters with the combined data. 
In this case, we perform the minimization procedure 
using Eq. (|4*3"|) without Xp r ior> the Xmin value is 73.50 
and (Xdof = 1-10). The best-fit value corresponds to the 
point (0.32, 0.34) in the fi m — n parameter space, and the 
confidence contours are shown in Fig[T2J We can see that 
the cosmological constant A lies in the 2-cr level. 

We collect our results for the several cases analyzed above 
in Tablc- HVl One can see that our results are consistent 
with the result obtained by Bento et al. [5^| , Gong and 
Duan [55|], Zhu et al. [57j]- However in our analysis the 
cosmological constant lies in between 1-cr and 2-cr. 



Cardassian's redshift 

With the knowlegde of these best-fit parameters we 
calculate z car d given by Eq. (jSJ) . Notice that this redshift 
is the one where the universe enters the phase of 
accelerate expansion driven by matter alone. This is 
the so-called Cardassian era. In general, it does not 
correspond to transition era to the accelerate expansion. 
Our results for the first three cases give 0.04 < z ca rd < 
0.56, in the most realistic case, where we considered 
the minimization procedure with the modified-;^ 2 for the 
GRBs data alone. For the combined data we obtain 
0.41 < z car d < 0.55. These results agree with the first 
initial values to z C ard given in the Table-I of the first 
Cardassian model proposed by Fresse and Lewis [l3l |. 
which was obtained by using observations of SNIa and 
the CMB. For each case, the values obtained for this 
redshift are given in Column (2) in the Table-TVl 



VI. 



OTHER OBSERVATIONAL ASPECTS 



Deceleration Parameter 



Finally, we shall analyze the influence of the best- 
fit parameters on various cosmological parameters, 
including the Cardassian redshift of transition, which is 
one of the main observational features of that model. 



This parameter was introduced to characterize the 
dynamical evolution of the universe based on the sign 
of the second time derivative of the scale factor. It is 
defined as 
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TABLE V: Other observational properties 



aa 
' ~7a 



(46) 



The negative sign was put by-hand to come along 
with ancient ideas and "observations" suggesting that 
the universe was slowing down its expansion since the 
Big Bang. As we are interested in expressing it as a 
function of both the Hubble parameter and cosmological 
redshift we can easily rewrite Eq. (|46p in the form: 



(47) 



l (l + z)dE 2 (z) 
q{z) = 2E*(z) dz - 1 



It is clear that a negative value of this parameter stands 
for an accelerating universe, and conversely, a positive 
value implies a decelerating universe. For the Cardassian 
model we obtain 



?(*) 



3(1: 

2 (1 



n)n(z) 



k{z)) 



(48) 



where n(z) = (j^j (1 + z)-^ 1 -^ . In Fig[13]we show 

the deceleration parameter as a function of the redshift 
for the values of the parameters obtained with the best-fit 
analysis. 

The value qo = q(z — 0) indicates the expansion rythm 
today. In this case we obtain 



1 3, 
q ° = 2 ~ 2^ ' x ' 



(49) 



The values obtained for this parameter lay in the range 
—0.91 < qo < —0.05 and are shown in Column (3) of 
Table-El 

The transition redshift z acc can be obtained by solving 
the equation q(z) — 0. Our result agrees with the SNIa 
observations 0.18 < z acc < 0.68. The corresponding 
values are shown in Column (4) of Table-TVl 



Age of universe 

Finally, we discuss the age-redshift relation in an 
alternatively form which is independent of the Hubble 
parameter today. To this purpose we use the 
measurement of the quantity Ho to'. 



Hoto 



dz 



(l + z)E(z) 



(50) 



where the dimensionless function E{z) is obtained from 
Eq.(fl2]). Our result shows that 0.86 < H t < 0.94 (see 
Column (5) in Table-TV]) . In order to verify our result 
we have compared with several other estimates of the 
same parameter. For instance, the CMB alone contrains 



Case c 


^cara 


fin 




rr , 

■£* ^0 


(1) 


(2) 


(3) 


(4) 


(5) 


exact Hq 


0.04 


-0.91 


0.37 


0.90 


exact Ho, prior Q m 


0.56 


-0.29 


0.68 


0.94 


Prior Ho 


0.19 


-0.20 


0.34 


0.88 


Prior Ho+Tl 


0.41 


-0.08 


0.24 


0.90 


Prior H + A 


0.55 


-0.05 


0.18 


0.86 


Prior H Q + TZ+ A 


0.46 


-0.18 


0.46 


0.89 



c Here we refer to a prior as a free parameter with Gaussian Prior 
distribution. 



to = 13, 7 ±0.25 Gyr[47j, while the Globular Cluster 
Age imposes to = 12.6^| 4 for the age of universe. 
If we also assume that Hq — /i(9.77813Gyr) , with 
h = 0.72±0.08 from the HST Key Project, and supposing 
that the Hq and to measurements are uncorrelated, we 
obtain for the product H n t the following range: 0.79 < 
H t < 1.27 and 0.67 < H to < 1.31, respectively. 
Clearly, one can see that our result is consistent with 
those estimates for the age of the universe. 



VII. DISCUSSION AND CONCLUSIONS 

In the Cardassian model, the equivalent of the dark 
energy density arises from a modification of Friedmann's 
equations through the introduction of a new function of 
the matter energy density in the form of an additional 
nonlinear term of mass. The outcoming model of the 
universe is then flat, matter dominated and accelerating. 

Meanwhile, in order to understand the expansion 
history of the universe one needs to be able to trace 
it back in time up to very high redshifts and construct 
its Hubble diagram from light-emitting cosmic sources 
laying along the pathway to the Big Bang. In this 
perspective, GRBs offer a mean to extend up to redshifts 
z > 6 the HD already built on SNIa observations 
z < 2. GRBs are now known to have several light- 
curves (in different energy bands) and various spectral 
properties from which the luminosity of each burst can be 
calculated, once it is calibrated for a specific cosmology. 
This last procedure has been proved to be self-consistent 
pol ]. and it turns GRBs useful standard candles for 
cosmographic studies. 

Following the same procedure used by Schaefer in 
BS2007, we made a simultaneous calibration with the 
Cardassian Cosmology of those GRB events that had 
their redshifts properly estimated. This allows us 
to construct the GRBs HD after using five empirical 
relations involving luminosity distance indicators of these 
GRBs. The results of the calibration procedure are 
shown in Figs. 2-6. The Fig.l states clearly that our 
method is self-consistent after comparing our results for 
the HD of GRBs calibrated with the concordance model 
with the one obtained by Schaefer 20]. The HD for 
the GRBs calibrated with the Cardassian model and 
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1 2 3 4 5 6 7 
Redshift (z) 

(a) Deceleration parameter for GRBs data alone 




Redshift (z) 



(b) Deceleration parameter for combined GRBs + 
CMB data 



FIG. 13: Deceleration parameters for Cardassian model as a function of the redshift for the values of best-fit 
solution. The horizontal line, labeled (d)/(a) (q(z) = 0), divide the region with (a) acelerating or (d) decelerating 
expansion at a given redshift. Clearly, we can see a transition from decelerating-to- accelerating universe in the range 
0.18 < z < 0.7, which is approximately the redshift range inferred from SNIa observations. 



its comparison with the A-CDM is presented in Fig. 7. 
We can see that the data of observed GRB events are 
consistent with both models. The small difference can 
be considered negligible if one takes into account that 
(a) the x 2 fo r both calculations are quite similar, and (b) 
the error bar of each GRB event is still quite large. 

We have also performed a detailed statistical analysis 
and constructed the confidence-level contours from these 
GRBs data after imposing several constraints on the 
cosmological parameters to be used for the analysis. 
First, considering exact knowledge of the Hubble 
parameter today, Hq = 72kms~ 1 Mpc _1 , the best- 
fit value obtained corresponds to the point (0.45, —0.7) 
in the f2 m — n parameter space. This point is shown 
in the associate confidence- level contour in Fig. 8. The 
same figure shows the case where one considers also a 
uniform prior knowledge on the matter density (f2 m = 
0.28 ±0.04), which corresponds to the point (0.28,0.27). 
In Fig. 9 we present the results of the analysis of the 
most specific case, that is, if one considers the lone 
GRBs dataset. We imposed a prior knowledge on the 
Hubble parameter today (iJ n = 72 ± 8 kms~ 1 Mpc _1 ) , 
and the best-fit value obtained corresponds to the point 
(0.40,0.26). This result also agrees with the first one 
obtained by Schaefer (2007). 

We have also analyzed the bounds from GRBs when 
combined with constraints provided by the CMB and 
BAO datasets. Considering the same prior knowledge on 
the Hubble parameter today, as we did before, the best- 
fit values obtained correspond to the points (0.35, 0.39) in 
the case of GRBs + K (Fig.10), (0.33, 0.45) in the case of 



GRBs + A (Fig. 11), and (0.32, 0.34) in the case of GRBs 
+1Z+ A (Fig. 12). The best-fit parameters obtained from 
this analysis are in agreement with those obtained from 
previous investigations that used SNIa data combined 
with CMB and BAO constraints. 

Besides, with the knowledge of the best-fit parameters, 
we analyze their influence on other observational 
properties of this model. For the Cardassian redshift, 
our results for the first three cases (with fixed value of 
H ) give 0.04 < z car( j < 0.56. For the combined data we 
obtained 0.41 < z cax d < 0.55. These results agree with 
the very first initial values of z C ard calculated by Freese 
and Lewis [l3j , which was obtained by using observations 
of SNIa and the CMB. For the deceleration parameter, 
the values obtained lay in the range —0.91 < qo < —0.05. 
For the transition redshift z acc we obtained 0.18 < 
^acc < 0.68, which agrees with that inferred from SNIa 
observations. Finally, we obtained a couple of bounds for 
the age-redshift product H t : (a) 0.79 < H to < 1-27 if 
one considers t = 13.7 ± 0.25 Gyr from CMB and 
(6) 0.67 < H t < 1.31, if one considers t = 12.6^ 4 
Gyr from the Globular Cluster Age Our results 

are consistent with those estimates for the age of the 
universe. 

In summary, the GRBs HD has been built to 
demonstrate the reliability of GRBs as an independent 
tool to check the consistency of most current cosmological 
models over distance scales not allowed to SNIa. This 
diagram for GRBs shows a great potential to make 
the GRBs observations a complementary tool to SNIa, 
large scale structure, baryon acoustic oscillations, an 
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perhaps to CMB observations. In the near future, new 
studies based on these novel relations promise a major 
breakthrough in using GRBs to do precision cosmology. 
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